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ABSTRACT 

We present aluminium abundances for a sample of about 100 red giant stars in each of the Galactic globular clusters 47 Tuc (NGC 104) and M 4 
(NGC 6121). We have derived homogeneous abundances from intermediate-resolution FLAMES/GIRAFFE spectra. Aluminium abundances 
are from the strong doublet Al i 8772-8773 A as in previous works done for giants in NGC 6752 and NGC 1851, and nitrogen abundances are 
extracted from a large number of features of the CN molecules, by assuming a suitable carbon abundance. We added previous homogeneous 
abundances of O and Na and newly derived abundances of Mg and Si for our samples of 83 stars in M 4 and 116 stars in 47 Tuc to obtain 
the full set of elements from proton-capture reactions produced by different stellar generations in these clusters. By simultaneously studying 
the Ne-Na and Mg-Al cycles of H-burning at high temperature our main aims are to understand the nature of the polluters at work in the first 
generation and to ascertain whether the second generation of cluster stars was formed in one or, rather, several episodes of star formation. Our 
data confirm that in M 4 only two stellar populations are visible. On the other hand, for 47 Tuc a cluster analysis performed on our full dataset 
suggests that at least three distinct groups of stars are present on the giant branch. The abundances of O, Na, Mg and Al in the intermediate 
group can be produced within a pollution scenario; results for N are ambiguous, depending on the C abundance we adopt for the three groups. 

Key words. Stars: abundances - Stars: atmospheres - Stars: Population II - Galaxy: globular clusters - Galaxy: globular clusters: individual: 
NGC 104 - Galaxy: globular clusters: individual: NGC 6121 
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1. Introduction 

The new paradigm for globular cluster (GCs) formation is cur- 
rently well assessed in its general aspects. Several decades of 
spectroscopic observations as well as years of high precision 
imaging mainly with HST have shown that a typical pattern 
is present within each galactic GC investigated so far (see the 
recent reviews by Gratton, Sneden and Carretta 2004; Martell 
2011; Gratton, Carretta and Bragaglia 2012a and references 
therein). A first generation of stars is found side by side with 
another, slightly younger stellar generation whose components 
are unambiguously formed by ejecta processed in the nuclear 
H-burning at high temperature within massive stars of the first 
generation in each cluster (see e.g. Gratton et al. 2001 for the 
first clear-cut observational constraint - followed by Ramirez 
and Cohen 2002, Carretta et al. 2004 -, and Ventura et al 2001 
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* Based on observations collected at ESO telescopes under pro- 
gramme 085.D-0205 and on public data from the ESO/ST-ECF 
Science Archive Facility. 

** Tables 2 and 3 are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+k/???/??? 



and Decressin et al 2007 for different views on the nature of 
these early polluters of the intra-cluster medium). 

The concept of GCs being good examples of Simple Stellar 
Populations (SSP) has been eventually replaced by their being 
made of Multiple Stellar Populations (MSP) whose composi- 
tions are marked by a typical pattern: light elements show anti- 
correlations and correlations among species produced (such as 
Na, Al, Si) and destroyed (such as O, Mg, F) in the network of 
proton-capture reactions responsible for the nuclear process- 
ing of the polluting matter (e.g. Denisenkov & Denisenkova 
1989; Langer et al. 1993, Arnould et al. 1999). In almost all 
GCs studied so far, ranging from the most massive one to the 
smallest objects in the Harris catalogue (1996, and following 
updates), this pattern is retrieved and it is so typical a signature 
of GCs that a new recent definition of bona fide globular has 
been proposed entirely based on the presence of the most no- 
table signature, the Na-O anti-correlation (Carretta et al. 2010), 
discovered by the Lick-Texas group in the 1990s (see the re- 
view by Kraft 1994) and extensively studied in many clusters 
(see the updated review by Gratton et al. 2012a and references 
therein.). 

While the first level of the game is clear, several details still 
are left in the poorly known age when GCs formed about a 
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Hubble time ago. Existent data are not enough to uncover yet 
the true nature of the polluters that provided the material for the 
building up of the second generations. The most favourite can- 
didates could be both intermediate-mass AGB stars (Ventura et 
al. 2001) or very massive, fast rotating stars (Decressin et al. 
2007), although alternative ones have been presented, such as 
massive binaries (de Mink et al. 2009). Moreover, the distribu- 
tions of stars along e.g. the Na-O anti-correlations do require 
that the matter composing the second stellar generation must 
be diluted in some amount with material that did not experi- 
ence hot nuclear burning. This rather strong constraint stems 
from both observations (the relations between proton-capture 
elements such as O and Na and the fragile element Li, easily 
destroyed at temperatures much lower than those involved in 
the hot H-burning, see Prantzos and Charbonnel 2006 and the 
discussion in Gratton et al. 2012a) and from theory (the yields 
from AGBs would result necessarily into a correlation, unless 
some amount of dilution is involved, D'Ercole et al. 201 1). 

Leaving aside for the moment several problems related to 
the nature of the diluting matter (D'Ercole et al. 2011), it is 
important to ascertain whether a model of dilution does really 
work and what are the observational consequences concerning 
the formation of multiple stellar generations in GCs. Two sce- 
narios can be envisioned: 

1 . After the formation of the first generation, a single mixture 
of matter from polluters diluted with unprocessed matter 
produces a second generation of stars distributed along the 
Na-O and Al-Mg anti-correlations. Since what is changing 
along these distributions is the degree of dilution, we ex- 
pect that elements destroyed in the proton-capture reactions 
(namely O and Mg) or species enhanced in these processes 
(i.e. Na and Al) show a direct, linear proportionality. 

2. On the other hand, several episodes of star formation might 
have occurred in the proto-GC after the formation of the 
first generation. Each one of them might have been char- 
acterized by a typical polluter with a given average mass. 
In this case, we expect that the elements are not necessar- 
ily produced or destroyed in linear proportion with each 
other. The distributions might still grossly resemble an anti- 
correlation, but different groups might be distinguished. 

Ideally, we would like to have for these studies a database 
composed of many stars with detailed abundances of several 
light elements. This requires extensive sets of very high reso- 
lution and S/N spectra with large spectral coverage, such as the 
FLAMES/UVES spectra obtained by Marino et al. (2008) who 
provided abundances for a good number of proton-capture ele- 
ments (O, Na, Mg, and Al for 93, 105, 105, and 87 stars, respec- 
tively) for M 4. When compared to our analysis for this cluster 
based on moderate resolution GIRAFFE spectra (Carretta et 
al. 2009a), their analysis show a clumpy distribution that was 
masked by noise on our lower quality data, despite the small 
star-to-star errors derived by our tailored procedure. Unluckily, 
it is not easy to obtain such high quality material for a large 
number of clusters because such observations are very expen- 
sive in telescope time. 

On the other hand, much insight can still be obtained by 
extending the spectral range of our moderate resolution survey 



to other spectral ranges, allowing to observe transitions for el- 
ements whose data were missing in our original spectra. This 
can be obtained with a very limited additional observing time 
exploiting the GIRAFFE spectrograph, which allows simulta- 
neous observation of more than a hundred star in a cluster, and 
can then be done easily for several clusters. Using this tech- 
nique, we showed e.g. that the globular cluster NGC 6752 is 
made of three discrete populations (Carretta et al. 2012a; note 
that the presence of three groups of stars in NGC 6752 is quite 
clear e.g. in the results obtained by Yong et al. 2008, but they 
did not make any comment about this in their paper). Moreover, 
we could show that the composition of the intermediate pop- 
ulation of that cluster cannot be reproduced by simply mixing 
suitable amounts of the material which made the two other pop- 
ulations. This clearly points towards the second scenario men- 
tioned above, possibly supporting the finding by Yong et al. 
(2008) that contributions from both AGB and massive stars are 
required to explain the correlations in NGC 6752. 

The present paper presents similar results on Al for two 
additional GCs (47 Tuc=NGC 104 and M 4=NGC 6121) for 
which we have already produced Na, O abundances in our 
FLAMES survey. The analysis is similar to that presented in 
Carretta et al. (2012a) for NGC 6752 and Carretta et al. (2012b) 
for NGC 1851. 

While for M 4 our data largely overlap with those of Marino 
et al. (2008), for 47 Tuc, this is the largest sample of RGB 
stars with homogeneous Na, O and Al abundances that we are 
aware of; the previous largest sample observed at medium-high 
resolution is represented by the 14 stars with UVES spectra 
analysed by our group in Carretta et al. (2009b), for Al, and 
Carretta et al. (2009a) concerning Na and O abundances. 



2. Observations and analysis 

The present work is based on two exposures, the first of 2700 
sec at airmass 1.487 (seeing=1.23") acquired on 2010, July 30 
for 47 Tuc, and the second of 1600 sec at airmass 1.240 (see- 
ing=1.79") taken on 2010, April 14 with FLAMES @VLT-UT2 
and the high resolution grating HR21. The resolution is 17,300 
at the center of the spectra and the spectral range goes from 
about 8484 to about 9001 A, including the Al i doublet at 8772- 
73 A, which is the main feature we were interested in. The po- 
sitioning of fibers used the same configurations employed to 
observe the strong doublet of Na (at 5682 and 5688 A) with the 
grating HR1 1 in Carretta et al. (2009a) to maximize the number 
of stars in the two clusters along the expected Na-Al correla- 
tion. We observed 116 giants in 47 Tuc in the magnitude range 

V = 12.25 - 14.57 and 83 RGB stars in the magnitude range 

V = 11.63 - 14.08 in M 4. Coordinates and magnitudes for 
all stars in our samples can be retrieved in the on-line tables of 
Carretta et al. (2009a). 

Data reduction was performed by the ESO personnel 
through the dedicated ESO FLAMES -GIRAFFE pipeline. 
The output bias-corrected, flat-fielded, one-dimensional and 
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wavelength-calibrated spectra were sky subtracted and shifted 
to zero radial velocity using 

The signal-to-noise (S/N) ratio of individual spectra is 
listed in Table 2 and Table 3 (only available on-line) for 47 Tuc 
and M 4, respectively. The median S/N is 249 for spectra of 
stars in 47 Tuc and 168 for giants observed in M 4. 

3. Analysis and derived abundances 

The atmospheric parameters required for the analysis were sim- 
ply adopted from Carretta et al. (2009a) where a full abun- 
dance analysis was done for all the stars targeted in the present 
study. This approach has two main advantages: first, it provides 
a full set of homogeneous parameters, with very small star- 
to-star errorfl useful to minimise down a possible smearing 
along the observed anti-correlations due only to errors asso- 
ciated to the analysis. Moreover, the number of lines of iron 
with well known atomic parameters is not high in the spectral 
range of the grating HR21, whereas the [Fe/Flj^| values derived 
in Carretta et al. (2009a) rest on several iron lines. Second, O 
and Na abundances, homogeneously derived using the same 
atmospheric parameters, can be safely coupled with the new 
abundances of Al, Mg, and Si. 

The complete list of atmospheric parameters (effective tem- 
perature T e ff, surface gravity, metallicity and microturbulent 
velocity), as well as the typical associated errors, can be re- 
trieved from the on-line tables in Carretta et al. (2009a) for all 
stars in 47 Tuc and M 4 analysed in the present study. 

3.1. Nitrogen 

CN lines are ubiquitous in the observed spectral range, and 
some of them could also contaminate the Al lines. Thus, the 
first step is to reproduce the CN features as much accurately 
as possible. The problem is more severe in metal-rich clusters 
such as 47 Tuc ([Fe/H]= -0.77 dex, Carretta et al. 2009c), but 
some degree of contamination should also be expected in the 
more metal-poor cluster M 4 ([Fe/H]= -1.17 dex, Carretta et 
al. 2009c). 

To deal with this aspect, we used an improved version of the 
procedure adopted in Carretta et al. (2012b). We first obtained 
a coadded spectrum with very high S/N value by summing up 
the ten spectra of 47 Tuc stars with the highest S/N. On this 
master spectrum, we individuated eight regions apparently free 
of lines to be used to derive a local reference continuum and 1 8 

1 IRAF is the Image Reduction and Analysis Facility, a general pur- 
pose software system for the reduction and analysis of astronomical 
data. IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in Tucson, 
Arizona. NOAO is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc. under cooperative agreement 
with the National Science Foundation. 

2 Typical internal errors are 4-6 K in temperature, 0.04 dex in grav- 
ity, 0.03 dex in [Fe/H] and 0.02-0.03 km/sec in microturbulent veloc- 
ity, see Table A.2 in Carretta et al. 2009a. 

3 We adopt the usual spectroscopic notation, i.e. for any given 
species X, [X]= log e(X) 8taI - log e(X) a and log e(X) = log (N x /N a ) + 
12.0 for absolute number density abundances. 



Table 1. List of regions used to estimate local continuum and 
fluxes of CN bands 
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CN 
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1.79 


cont 
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cont 


8791.26 


8791.46 


0.20 


cont 
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8792.12 


0.38 


cont 


8794.76 


8795.32 


0.56 


cont 


8834.69 


8835.14 


0.45 



spectral regions dominated by CN bands. Definitions of these 
spectral regions and of those for the local reference continuum 
are given in Tab. Q] 

In the second step, synthetic spectra were computed us- 
ing the Kurucz (1993) grid of model atmospheres (with the 
overshooting option switched on) and the same line list used 
in Carretta et al. (2012b), using the atmospheric parameters 
suitable for each individual star in 47 Tuc. Three spectra were 
computed assuming [C/Fe]=-0.6 dex. and three evenly spaced 
values of nitrogen, [N/Fe]=-0.5, +0.25 and +1.0 dex, with the 
package ROSA (Gratton 1988). To each of the 18 CN features 
in each star we applied the procedure explained in Carretta et 
al. (2012b): the average flux within the in-line region of the 
feature was measured and a weighted reference continuum was 
derived from all the 8 off-line regions using weights equal to the 
width of each region. Whenever a defect in the spectrum was 
found on a continuum region from visual inspection, that re- 
gion was excluded from the computation of the reference local 
continuum. The N abundance from each feature was then ob- 
tained by comparing the normalized flux with the fluxes mea- 
sured in the same way on the three synthetic spectra. 

The above procedure was repeated for each of the 116 stars 
observed in 47 Tuc and for each star an average [N/Fe] value 
was derived. The spectral regions of each CN feature were 
inspected visually and the relative feature discarded from the 
average whenever a spike was found to fall on the feature. 
The final abundances of N were obtained after applying a kcr- 
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Fig. 1. Left: comparison of our derived [N/Fe] ratios with the 
values of the DDO index from Norris and Freeman (1979) 
in 47 Tuc. Right: comparison of our N abundances with the 
S(3839) index for the strength of the ,13883 CN band from the 
compilation by Smith and Briley (2005) in M 4. 



clipping at 2.5 cr to the average abundance from individual 
features in each star. These final abundances in 47 Tuc rest 
on average on 17 features. We found a mean value [N/Fe] = 
+ 1 .038 + 0.01 1 dex (rms=0.1 19 dex, 1 16 stars). 

This approach was adopted also for the 83 stars observed in 
M 4; also in this case, the mean N abundance was derived from 
an average of 17 CN features, finding [N/Fe]= +1.052 ± 0.017 
dex (rms=0.152 dex, 83 stars). The larger scatter observed in 
M 4 is mainly due to the giants in this more metal-poor cluster 
being on average 200 K warmer than giants in 47 Tuc (and only 
minimally due to the lower S/N of the spectra). The molecular 
features are more difficult to measure and the warmest stars are 
thus responsible for most of the observed scatter. 

Individual values of N abundances are listed in Tab. 2 and 
Tab. 3 for 47 Tuc and M 4, respectively. We recall that these 
[N/Fe] were obtained assuming [C/Fe]=-0.6. Were the C abun- 
dance for a given star be different from this value, [N/Fe] 
should be modified in such a way that the sum [C/Fe] + [N/Fe] 
remains constant, because the concentration of CN depends on 
the product of the C and N concentrations. 

The C abundances are expected to change as a function 
of temperature due to deep mixing. We estimated from C and 
N abundances measured in red giants of M 4 by Ivans et al. 
(1999) that over the temperature range of our sample the maxi- 
mum variation in C expected as a function of the luminosity is 
about 0.25 dex. The effect of star-to-star variations in C abun- 
dances at a given luminosity, related to the multiple population 
phenomenon, will be discussed later. The impact is likely very 
important in the case of 47 Tuc. For M 4, it seems less of a 
problem. For instance, we have sixteen stars in common with 
Smith & Briley (2005) that obtained C abundances from low 
resolution spectra. If we divide these stars in 6 Na-poor and 10 
Na-rich stars (at [Na/Fe]=0.4), we find [C/Fe]=-0.57 + 0.04 
(r.m.s.=0.10) and [C/Fe]=-0.64 + 0.06 dex (r.m.s=0.20 dex) 
for the two groups, both compatible with our assumption of 
[C/Fe]=-0.6. 



We tested the adopted procedure on the spectrum of 
Arcturufl We obtained [N/Fe] = 0.68 adopting [C/Fe]=0 and 
the atmospheric parameters from Hinkle et al. (2000) and 
Carlberg et al. (2012): T eff = 4300 K, logg = 1.70, [Fe/H]=- 
0.55, and v, = 1 .64). This result is very similar to that obtained 
by D'Orazi et al. (201 1) with very similar atmospheric param- 
eters. 

For a check, we searched the literature for analysis with 
published values of N or CN (or related indexes) in these two 
clusters. We found 1 1 RGB stars of 47 Tuc in common with 
the sample studied by Norris and Freeman (1979), and we 
measured N abundances for 9 of them; we also have 29 stars 
in common with the compilation of CN values by Smith and 
Briley (2005) in M 4, with our N values available for 18 of 
them. In Fig.Q]we show the comparison of our derived N abun- 
dances with the values of the DDO index C(4142) in 47 Tuc 
(left panel) and of the re-homogenized index S(3839) by Smith 
and Briley (2005) in M 4 (right panel). In both cases, with our 
assumptions about the carbon abundance, our values are quite 
well correlated with these literature values. We conclude that 
our values provide a good estimate of the N content in the pro- 
gram stars. 



3.2. Aluminium, magnesium and silicon 

A similar procedure was used to derive Al abundances (see also 
Carretta et al. 2012b). We used for each star three synthetic 
spectra computed with the appropriate atmospheric parame- 
ters (from Carretta et al. 2009a) and the N abundances derived 
as described above. The fluxes for determining Al were mea- 
sured in the in-line region from 8772.4 to 8774.7 A while two 
slightly different regions were used in the two clusters to derive 
the local reference continuum owing to the different metallic- 
ity of these GCs. The selected regions were 8771.28-8771.70 A 
and 8776.54-8777.46 A in 47 Tuc, and 8768.90-8770.00 A and 
8775.20-8776.10 A in M 4. Again, we obtained a final contin- 
uum by averaging the contribution from the two regions with 
weights given by the number of pixels (i.e. the widths of these 
intervals). 

Abundances of Al were then derived by interpolating the 
normalized fluxes in the in-line Al region among those obtained 
from the synthetic spectra, computed for [Al/Fe]=0.0, 0.75, and 
1 .5 dex. As above, should one of the reference local continuum 
regions or one of the Al lines be affected by a defect, following 
visual inspection, the contaminated region would be dropped 
from the computation. 

Abundances of Mg and Si were obtained from the same set 
of lines described in Carretta et al. (2012a), using atmospheric 
parameters derived in Carretta et al. (2009a). 

Final Al, Mg, and Si abundances for giants in 47 Tuc and 
M 4 are listed in Tab. 2 and Tab.3, respectively. 



4 We retrieved this spectrum in the relevant wavelength 
range from the ESO Advanced Data Product Archive 
ihttp : // archive.eso.org/eso/eso.archive_adp.html). 
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3.3. Errors 

Star-to-star errors in the adopted atmospheric parameters were 
estimated in Carretta et al. (2009a) and are shown in Tab. [4] 
To translate these errors into errors on the Al abundances, we 
need the sensitivities of Al to changes in the adopted parame- 
ters. These were obtained by repeating the analysis for the stars 
29490 (T eff = 4584 K) in 47 Tuc and star 30450 (T eff = 4622 
K) in M 41 and changing a parameter each time, holding the 
others fixed at their original value. The amount of the varia- 
tions and the sensitivity of each parameter to these changes are 
listed in Tab. |4] 

The star-to-star (internal) errors in the Al abundances are 
obtained by summing in quadrature all the contributions, in- 
cluding the impact of errors in the derivation of N abundances 
from CN features. 

We attached an error due to flux measurement to the Al 
abundances. This was evaluated, as in Carretta et al. (2012b) 
estimating photometric errors from the S/N ratio of the spectra 
of stars and from the width within each of the reference contin- 
uum and in-line regions (and then the number of independent 
pixels used). The error in Al abundances is then obtained by 
comparing these Al abundances with those derived considering 
a new value of the Al line strength index that is the sum of the 
original values and its error. The errors from these procedure 
are listed for each star in Tab.2 and Tab.3; on average, they are 
+0.013+0.001 dex from 1 16 stars in 47 Tuc and +0.019+0.001 
dex from 83 giants in M 4. 

Moreover, we applied the same approach to one of the CN 
features used to estimate the N abundance. The average error 
is 0.012 dex in 47 Tuc and +0.028 dex in M 4. These values 
were assumed as conservative estimates in the error budget, 
since final N abundances were obtained from a large number 
of features. 

Summing up all the above contributions the typical internal 
errors in the derived Al abundances are 0.042 dex and 0.036 
dex for 47 Tuc and M 4, respectively. Typical internal errors in 
[Mg/Fe] and [Si/Fe], estimated according to the usual proce- 
dure (Carretta et al. 2009a), are 0.10 dex and 0.06 dex, respec- 
tively, for both clusters. 

4. Results 

The pattern of abundances of proton-capture elements among 
red giants in 47 Tuc is summarized in Fig.|2]and Fig. [3] 

The Na-O anti-correlation in 47 Tuc was already known to 
be rather short, when considering its large total mass, which 
is one of the main parameters driving the extent of the anti- 
correlation (Carretta et al. 2010). The new set of light elements 
added in the present work confirms this trend: Mg and Si do 
not show large star-to-star variations, the average values be- 
ing [Mg/Fe] = +0.532 ± 0.007 dex (rms=0.079, 147 stars) and 
[Si/Fe] = +0.440 ± 0.005 dex (rms=0.065, 147 stars JB The alu- 



5 These stars were chosen because they have an effective tempera- 
ture in the middle of the temperature range of their respective samples. 

6 The number of stars where Mg and Si is measured is higher than 
those with Al abundances since transitions of Mg and Si are available 
in the spectral range of both the gratings HR1 1 and HR13, whereas Al 
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Fig. 2. Run of the abundance ratios of proton-capture elements 
O, Na, Mg, Al and Si among red giants in 47 Tuc, from Carretta 
et al. (2009a) and the present study. Upper limits in O are indi- 
cated by arrows. Typical star-to-star errors are also indicated. 
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Fig. 3. Run of the nitrogen abundances as a function of O, Na, 
Mg, and Al for our sample of RGB stars in 47 Tuc. 
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Table 4. Sensitivities of Al to variations in the atmospheric parameters and to errors in fluxes, and errors in abundances [A/Fe] 
for stars in 47 Tuc and M 4 



47 Tuc 



Element 


[N/Fe] 


Tell' 


l°g# 


[A/H] 


v, 


flux 


Total 






(K) 


(dex) 


(dex) 


kms -1 


(dex) 


Internal 


Variation 


+0.200 


50 


0.200 


0.100 


0.10 






Internal 


+0.012 


6 


0.020 


0.032 


0.11 


0.013 





[Al/Fe]i 0.077 -0.032 0.000 0.115 +0.013 0.042 



M4 



Element 


[N/Fe] 


Tell' 


logs 


[A/H] 


v, 


flux 


Total 






(K) 


(dex) 


(dex) 


kms~' 


(dex) 


Internal 


Variation 


+0.200 


50 


0.200 


0.100 


0.10 






Internal 


+0.028 


4 


0.041 


0.025 


0.12 


0.019 




[Al/Fe]i 


0.013 


-0.031 


0.030 


0.104 


+0.012 




0.036 



minium abundances span a range of about 0.5 dex among giants 
in this cluster, with an average value of [Al/Fe] = +0.529+0.012 
dex (rms=0. 132, 1 16 stars). Clear trends of Al ratios correlated 
to Na and anti-correlated to O abundances show the presence 
in second generation stars of ejecta processed through the high 
temperature Mg-Al cycle. This finding is also confirmed by the 
correlation (shallow, but still statistically robust) between Si 
and Al, that by itself points toward temperature of H-burning 
in excess of ~ 65 x 10 6 K, where the leakage from the Mg- 
Al cycle on 28 Si is onset (see Yong et al. 2005, Carretta et al. 
2009b). 

As expected, the abundances of N derived from the trans- 
formation of O in the CNO cycle, are nicely correlated with el- 
ements enhanced in proton-capture reactions (such as Na, Al) 
and anti-correlated with oxygen (Fig. [3}. The apparent lack 
of a N-Mg anti-correlation is not a big source of concern, 
since the Mg variations are small. Our survey of Al and Mg 
abundances in several GCs, based on limited samples of RGB 
stars observed with FLAMES -UVES (Carretta et al. 2009b) 
show that Al-rich and Mg-depleted stars are present only in 
massive (NGC 2808, NGC 6388, NGC 6441) or metal-poor 
(NGC 6752) clusters, or both (NGC 7078=M 15). In such clus- 
ters a clear Mg-Al anticorrelation is observed even among main 
sequence stars (see e.g. Bragaglia et al. 2010 for NGC 2808). 
In GCs such as those under scrutiny, small and/or metal-rich, 
large star-to-star variations of Mg abundances are not expected, 
as shown by Carretta et al. (2009b) in a survey of limited sam- 
ples of giants in 15 GCs and by Marino et al. (2008) over about 
100 giants inM4. 

The analogous run of light elements in the globular cluster 
M 4 are shown in Fig. [4] and Fig. [5] In this small mass cluster 
the Na-O anti-correlation is found to be of moderate extension 
(Carretta et al. 2009a, Marino et al. 2008, Ivans et al. 1999) 
and also star-to-star variations in Mg and Si are not large. The 
average values from our analysis are [Mg/Fe] = +0.541 ± 0.005 
dex (rms=0.048, 103 stars) and [Si/Fe]= +0.540 + 0.003 dex 
(rms=0.034, 103 stars). As already indicated by previous anal- 
ysis (Marino et al. 2008, Carretta et al. 2009b), star-to-star vari- 
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from GIRAFFE spectra is available for only the subsample observed 
with the grating HR21. 



Fig. 4. The same as in Fig. |2]for M 4. 



ation in Al abundances are not large along the RGB in M 4. A 
clear Na-Al correlation is evident in our data, and in the UVES 
data of the large sample in Marino et al. (2008), whereas it can- 
not be seen in the limited sample of stars analyzed in Carretta 
et al. 2009b). On the other hand, no correlation between Si and 
Al can be appreciated in the present study. 

Concerning N abundances, we found a somewhat larger 
scatter with respect to the results of 47 Tuc. As explained 
above, this is largely due to the increased difficulty to measure 
CN features in the warmest stars of our sample in this more 
metal poor cluster. This is shown in Fig. [6] where [N/Fe] ratios 
are plotted as a function of the effective temperature: the scatter 
noticeably increases for stars at the warmest end of the sample, 
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Fig. 6. Abundances of N as a function of effective temperature 
in giants of M 4 



which are those accounting for most of the dispersion seen in 
Fig.H 

4.1. Comparison with previous studies 

We found in our sample 85 stars in common with the analysis 
of M 4 by Marino et al. (2008); the comparison between our 
abundances of Na, Al, Mg and Si is summarized in Fig. [7] 

Taking into account differences in the scale of atmospheric 
parameters (based on photometry in our case; eventually ad- 
justed using spectroscopic parameters, in the analysis made 
by M08), abundance indicators (M08 used the classical opti- 
cal doublet of Al i at 6696-98 A instead of the stronger doublet 
at 8772-73 A as in the present work), atomic parameters and 




0.2 0.4 0.6 0.8 
[Al/FcJ MOB 



0.2 OA 0.6 0.8 1 

[Si/FcJ M08 




0.2 0.4 0.6 
[Mg/Fe] M08 



0.5 0.5 1 
[Na/Fe] M08 



Fig. 7. Comparison of abundance ratios of Na, Al, Mg and Si 
for stars in our sample in common with the study by Marino et 
al. (2008). The solid line is the line of equality. 



the different resolution of the spectra (higher resolution UVES 
spectra for M08), the agreement is satisfactory. 

Concerning 47 Tuc, recently Gratton et al. (2012b) stud- 
ied a sample of 1 10 stars on the red horizontal branch (HB) in 
47 Tuc. The agreement in O and Na abundances is very good: 
the Na-O anti-correlation as observed along the RGB and the 
HB is almost indistinguishable (see Figure 5 in Gratton et al.), 
despite the different abundance indicators adopted in the two 
studies. However, for N we found that an offset exists, the bulk 
of [N/Fe] values being found between ~ 1.2 and ~ 2 dex for 
HB stars in 47 Tuc. Finally, some offset is also seen regarding 
Al abundances (Figure 14 in Gratton et al.), which are derived 
from different doublets (7835-36 A for HB stars and 8772-83 A 
here). 

4.2. Cluster analysis 

The large dataset (both in number of elements and in number 
of stars) available to us allows to perform a cluster analysis to 
find out whether stars observed along the relations among light 
elements do form discrete groups. As in Carretta et al. (2012a) 
this was done using the A:-means algorithm (Steinhaus 1956; 
MacQueen 1967) as implemented in the R statistical package 
(R Development Core Team 2011; http://www.R-project.org). 
For this exercise we selected the abundance ratios [Al/Fe], 
[O/Fe], [Na/Fe], [Mg/Fe], and [N/Fe], considering only stars 
in 47 Tuc and M 4 with all these elements simultaneosly deter- 
mined: 79 giants in 47 Tuc and 67 in M 4. The parameters were 
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1.5 2 2.5 3 
[(N + Na+Al)/Fe] 



Fig. 8. Upper panel: results of the cluster analysis on stars in 
47 Tuc with a full set of abundances of N, O, Na, Mg, Al. The 
input number of groups to be found is two. Lower panel: the 
same, but using three as the number of groups in input. 



weighted according to their internal errors, which are largest 
for O abundances, usually derived by only one line (the asso- 
ciated internal errors in O are for example about twice those 
relative to Na, see the Appendix in Carretta et al. 2009a). 

The number of subgroups to be found is a parameter that 
must be provided as input to the cluster analysis when using the 
£-means algorithm. The simplest hypothesis we could make is 
that two groups, i.e. stellar generations exist in the cluster. The 
first one has a primordial, typical supernovae-enriched abun- 
dance pattern; the second has a chemical composition modified 
by ejecta of a fraction of stars of the first generation (the pol- 
luters). 

The results are shown in the upper panel of Fig. [8] for stars 
in 47 Tuc; we used as x-axis a parameter including all the el- 
ements that are enhanced in the proton-capture reactions. This 
parameter quantifies how much the second generation is ex- 
treme. In ordinate we display the ratio [O/Fe] since it is the 
element most susceptible to depletion through these reactions, 
while Mg variations are, for example, usually very tiny. The av- 
erage abundances for each group are listed in Tab. [5] The algo- 
rithm basically finds in 47 Tuc two groups that corresponds to 
the two generations, assigning about half of stars to each group. 
A fraction of 50% of stars of the first generation is somewhat 
higher than the fraction of ~ 30% according to the definition of 
primordial P fraction given in Carretta et al. (2009a). 

However, looking at the upper panel in Fig. [8] we note 
that beside the rather clear separation between stars in the 
first (open circles) and in the second groups (filled circles) 
there is some evidence for a discontinuity in the distribu- 




1.5 2 2.5 3 
[(N + Na+Al) /Fe] 



Fig. 9. The same as in Fig. [8]for M 4. 



tion also within the first group, located approximatively at 
[(N+Na+Al)/Fe]=1.8 dex. To verify how sound is this feature, 
we ran again the cluster analysis, forcing the algorithm to find 
three groups. 

The results are presented in the lower panel of Fig. [8] the 
stars of the previous first generation are now split into two 
groups and the putative intermediate population detected above 
nicely falls into a unique group (filled triangles), apart from a 
single outlier. With this separation, the fraction of stars with 
primordial composition is now 24%, while the intermediate 
and extreme populations summed together reach 76% of the 
sample, reassuringly close to the typical fractions derived in 
Carretta et al. (2009a) using only O and Na abundances. This 
is not surprising, since the positioning of the GIRAFFE fibers 
with the HR21 grating in the present study, was done using 
the same configuration used for Na (HR11) by Carretta et al. 
(2009a). As explained in that paper, these configurations, due 
to the limitations of the Oz-Poz positioner of FLAMES, are not 
the best suited to study the radial dependence of the fraction of 
stars in the P, I and E components. This also makes difficult a 
direct comparison with methods based on photometric parame- 
ters, like those in Milone et al. (2012), more able to study radial 
distributions, but with different sensitivities on chemistry and 
age. For instance, we note that in Milone et al., three popula- 
tions were found among subgiants, but not among RGB stars. 
Also, we are not sure that the three populations found along 
the RGB in the present study by means of chemical composi- 
tion are the same found along the SGB. 

For M 4 the case is much simpler. The cluster analysis with 
only two groups as input separates quite well the first and sec- 
ond generation stars, as evident in the upper panel of Fig. [9] 
where the gap along the "secondness" parameter is very clear. 
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Table 5. Average values of abundance ratios of the different groups resulting from the cluster analysis in 47 Tuc and M 4 



47 Tuc: 2 groups 



group 




nr 


[N/Fe] 


[O/Fe] 


[Na/Fe] 


[Mg/Fel 


[Al/Fe] 


I 




41 


Tu.yo in u.Ui 


+n ?<s + o ni 


+o % + n n? 


4-0 S? + 01 

TU.Jz X U.Ul 




2 




38 


+ 1 11+001 


+0 00 + 0? 


+0 SQ + 0? 

TU.J7 LT_ VJ.VJZ, 


+0 S 2 + 01 


+0 <S4 + 01 


Al Tup 


— = 

: 3 groups 














1 
1 




1 Q 


+U.0 / ± U.U1 


_i_n ^ 1 _i_ n no 
+U..51 ± U.Uz 


, n oo i n no 
4-U.zo ± U.UZ 


+U. J 1 ± U.UZ 


, n /in , n ni 
+U.4U ± U.Ul 


2 




22 




+0 22 + 02 


4-0 42 4- 01 


4-0 S2 4- 02 


+ n 4Q + a 01 


3 




38 


+ 1.11 ±0.02 


+0.00 ± 0.01 


+0.59 ± 0.02 


+0.52 ±0.01 


+0.64 ± 0.01 


M4: 2 


groups 














group 




nr 


[N/Fe] 


[O/Fe] 


[Na/Fe] 


[Mg/Fe] 


[Al/Fe] 


1 




13 


+0.88 ± 0.01 


+0.35 ± 0.02 


+0.10 ±0.02 


+0.55 ±0.01 


+0.58 ± 0.03 


2 




54 


+ 1.12 ±0.03 


+0.19 ±0.03 


+0.45 ± 0.03 


+0.54 ±0.01 


+0.77 ± 0.02 


M4: 3 


groups 














1 




13 


+0.88 ±0.01 


+0.35 ± 0.02 


+0.10 ±0.02 


+0.55 ±0.01 


+0.58 ± 0.03 


2 




24 


+ 1.21 ±0.01 


+0.26 ± 0.02 


+0.42 ± 0.02 


+0.53 ±0.01 


+0.73 ± 0.03 


3 




30 


+ 1.05 ±0.02 


+0.13 ±0.03 


+0.47 ± 0.03 


+0.55 ±0.01 


+0.80 ± 0.02 



Hence, it is not surprising that the cluster analysis fails to find 
three distinct subgroups (lower panel in Fig.|9]i: in this case an 
"intermediate" group is found by the algorithm, but it is barely 
separated from the bulk of second generation stars. A look at 
the abundances listed in Tab. [5] shows that this putative inter- 
mediate population has N abundances too high, with respect to 
the extreme groups, and the likely explanation is that these are 
the warmest stars, responsible for the "flaring" of the scatter in 
N seen above. This is simply an effect of the analysis and not a 
true feature in this cluster. 

5. Discussion and conclusions 

The present study highlights the advantage of having data for a 
large set of elements involved in the network of proton-capture 
reactions of H-burning at high temperature. It is very useful 
to have a simultaneous look at all relations among light ele- 
ments because the depletions observed for O and Mg, as well as 
the enhancements detected in N, Na, and Al allow probing the 
stratigraphy within the H-burning region, reaching more and 
more hot layers. In turn, this translates into evidence of more 
massive stars being responsible as the source of intra-cluster 
pollution. 

This large harvest of data can then be exploited to under- 
stand whether the formation of globular clusters was charac- 
terized by a single second episode or multiple episodes of star 
formation. Of course, when the observational evidence points 
toward only two distinct groups, as found in M 4 in the present 
and previous analyses (Ivans et al. 1999, Marino et al. 2008, 
Carretta et al. 2009a), the question must necessarily remain 
unanswered and the scenario is undetermined. 

More interesting is the case when more than two groups are 
found. Leaving aside the peculiar case of u> Cen, the most no- 
ticeable case among normal GCs is NGC 2808. We know that 
this cluster shows three distinct main sequences (D'Antona et 
al. 2005; Piotto et al. 2007), has an unusual HB, bimodal and 
with a long blue tail (Walker 1999) that has been connected 



to different He enrichments (D'Antona and Caloi 2004), and 
shows a very extended Na-0 anti-correlation, with some evi- 
dence of multiple "blobs" (Carretta et al. 2006). 

Another more recent case is NGC 6752, where with our 
survey of Al abundances we found the evidence of three dis- 
tinct stellar populations with different chemical composition. 
Moreover, from the correlation between depleted species such 
as O and Mg (clear but not linear) we were able to show that the 
composition of the intermediate population cannot be obtained 
by simply mixing primordial matter with material having the 
chemical composition of the extreme population. 

The same approach can be followed to try to better un- 
derstand our present finding in 47 Tuc. Using a simple dilu- 
tion model (see Carretta et al. 2009a, Gratton et al. 2010) we 
searched for the dilution factors able to reproduce the interme- 
diate group by mixing the typical composition of stars from 
the other two groups. If the model is correct we should obtain 
the same factor, regardless of our choice of the element to be 
used. We derived that dilution factors 0.71, 0.64, and 0.63 (ba- 
sically, the same factor, within the uncertainties) are required 
for Al, O and Na, respectively. The Mg abundance variations 
are so small that virtually any factor applies. Our nominal N 
abundances, which were obtained assuming a constant value of 
[C/Fe] =-0.6 for all stars, require a dilution factor of 0.29 that 
is quite different from the values obtained for the other ele- 
ments. However, this result depends on the adoption of a con- 
stant [C/Fe] and a different conclusion is obtained by assuming 
that the C abundance is different in the three groups. This is 
not at all an ad hoc assumption, since several studies showed 
that the [C/Fe] ratio is higher by ~ 0.3 dex in CN-weak (Na- 
poor, O-rich) stars than in CN-strong (Na-rich, O-rich) stars 
(see Norris & Freeman 1982, Norris & Cottrell 1979; Smith et 
al. 1989; Brown et al. 1990; Briley 1997; Carretta et al. 2005). 
Were the C abundance in group 2 be also intermediate between 
those of groups 1 and 2, it would be possible to reproduce the 
strength of the CN lines with a N abundance compatible with 
a pollution scenario. Direct determination of C abundances for 
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the stars in our sample would solve this ambiguity, and clarify 
if dilution may explain the pattern seen in 47 Tuc. The same 
observation could also be used to find if the sum of C+N+O 
abundances is constant or not in 47 Tuc. We only note here 
that if we combine the C abundances for CN-weak and CN- 
strong RGB stars by Brown et al. (1990) and Briley (1997) 

([C/Fe] 0.1 and [C/Fe]~ -0.3, respectively) with our CN 

and O abundances, the sum of C+N+O abundances would be 
the same for the two extreme groups within 0. 1 dex. This re- 
sult agrees quite well with the study of the C+N+O sum in 
unevolved cluster stars in NGC 6752, NGC 6397, and 47 Tuc 
made by Carretta et al. (2005). 

Another problem in 47 Tuc is to find the progeny of the 
subgiant (SGB) stars C+N+O enhanced, the so called "SGII", 
hypothesised by Di Criscienzo et al. (2010) that are 10% of the 
stars. More investigations are needed to have a self-consistent 
picture on this cluster. 

A promising approach could be to look for the full set of 
abundances, from C, N, O to the light and heavy, neutron- 
capture elements, especially if measured in the same stars. 
However this is generally not available at the moment, so no 
stringent conclusions can be drawn. For instance, Worley et 
al. (2010) present a list of studies regarding heavy elements in 
47 Tuc. Only in a few of them, and only for a limited number 
of stars, abundances for the whole set of elements is available 
(e.g, 12 TO and SGB stars in James et al. 2004, when com- 
bined with Carretta et al. 2005; four bright giants in Brown & 
Wallerstein 1992). An interesting approach based on medium 
resolution spectra of a large number of targets is presented in 
Worley & Cottrel (2012), who obtained a measure of the CN 
excess for about 100 giants, with almost no overlap with our 
present sample. However, they did not measure O abundances 
and published results of the abundance analysis of other el- 
ements only for a small subsample of 13 bright stars. In all 
the above-mentioned studies there seems to be an homoge- 
neous pattern among elements produced by neutron-capture 
processes, and this conclusion seems to be a severe drawback 
for the hypothesis of a significant contribution of polluted mat- 
ter by low-mass AGB stars in this cluster. We stress however 
that all the existing samples are too small to permit a separa- 
tion of groups as we do in the present paper using the cluster 
analysis technique. Observational efforts should be devoted to 
acquiring the most complete set of key elements abundances in 
significant samples of stars. 

Our results show that 47 Tuc is the third normal globular 
cluster after NGC 2808 (D'Antona et al. 2005, Piotto et al. 
2007, Carretta et al. 2006) and NGC 6752 (Yong et al. 2008, 
Carretta et al. 2012a), where at least three groups of stars with 
homogeneous chemical composition (within each group) can 
be separated. The implication seems to be that for a fraction 
of GCs the scenario of cluster formation must include multi- 
ple bursts where second generation stars are produced. In each, 
the interplay of matter processed by polluters of different (ob- 
viously decreasing with time) average masses, and pristine ma- 
terial does combine to give the observed chemical pattern in 
each group. By studying the detailed abundance distribution of 
proton-capture elements we can then provide important con- 
straints to the theoretical models of cluster formation. 



All these ingredients must be carefully considered 
and tuned into existing and future theoretical models. 
Homogeneous analysis of several proton-capture elements in 
large sample of stars such as the one presented in this study 
may provide a fruitful way to contribute useful constraints to 
these models. 
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